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Abstract 



I A starquake mechanism for pulsar glitches is developed in the solid quark star 

j- ' model. It is found that the general glitch natures (i.e., the glitch amplitudes and the 

. time intervals) could be reproduced if solid quark matter, with high baryon density 

I but low temperature, has properties of shear modulus /i = 10^''~^^ erg/cm^ and 

' critical stress Uc = 10^^~^^ erg/cm^. The post-glitch behavior may represent a kind 

'nI" . of damped oscillations. 
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1 Introduction 



Pulsars are unique objects, with which all types of elemental interaction could 
be tested extremely. However, the most elementary question relevant is still 
open: What is the nature of pulsars? It is conventionally thought that pulsars 
are simply a kind of boring big "nuclei" — neutron starSjbut more and more 
attention is paid to the quark star model for pulsars (|l|; |2) since no convinc- 
ing work, neither in theory from first principles nor in observation, has con- 
firmed Baade-Zwicky's original idea that supernovae produce neutron stars. 
The bare quark surface is suggested to be a new probe for identifying quark 
stars with strangeness, and possible observational evidence for bare strange 
stars appears: the drifting subpulses of radio pulsars, ultra-high luminosity of 
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soft 7-ray repeaters, non-atomic thermal spectra of isolated "neutron" stars 
(j^. However, can the bare strange star model reproduce most of the general 
features of pulsars (especially glitches)? 

The observation of free precession in PSR B1828-11 (Q) and PSR B1642-03 (0) 
challenges astrophysicists today to re-consider the internal structure of radio 
pulsars The current model for glitches involves neutron superfluid vertex 
pinning and the consequent fluid dynamics. However, the pinning should be 
much weaker than predicted in the glitch models, otherwise the vortex pin- 
ning will damp out the precession on timescales being much smaller than ob- 
served. In addition the picture, that a neutron star core containing coexisting 
neutron vertices and proton flux tubes, is also inconsistent with observations 
of freely precessing pulsars (0). It is then supposed that the hydrodynamic 
forces presented in a precessing star are probably sufficient to unpin all of the 
vortices of the inner crust (8) since a definitive conclusion on the nature of 
vertex pinning has not been reached yet due to various uncertainties in the 
microscopic physics. But recently. Levin & D'Angelo (0) studied the magneto- 
hydrodynamic (MHD) coupling between the crust and the core of a rotating 
neutron star, and found that the precession of PSR B1828-11 should decay 
over a human lifetime. This well-defined MHD dissipation should certainly be 
important in order to test the stellar models. 



An alternative way to understand both glitch and free-precession could be 
through the suggestion that radio pulsars are solid quark stars (13 U). A solid 
quark star is just a rigid-like body, no damping occurs, and the sohd pulsar 
model may survive future observational tests if the free precession keeps the 
same over several tens of years. A neutron star could not be in a solid state, 
whereas a cold quark star could be. Such a solid state of quark matter could 
be very probably Skyrme-like 2 (0 13), the study of which may help us to 
understand dense quark matter with low temperature. 



Fluid strange-star (even with possible crusts) models were noted to be incon- 
sistent with the observations of pulsar glitches more than one decade ago (.15) . 
Modifications with the inclusion of possible stable particles to form a differen- 
tiated structure of so-called strange pulsars was also sug gest ed (IigI ). but is not 
popular because of a disbelief in the employed physics ()l7l : 0). However, can 
a fully solidified quark star proposed |lOl) really reproduce the glitch behav- 
iors observed? One negative issue is that giant glitches are generally not able 
to occur at an observed rate in a solid neutron star (18). Nonetheless, more 
strain energy could be stored in a solid quark star due to an almost uniform 
distribution of density (the density near the surface of a bare strange star is 
~ 4 X 10^^ g/cm^) and high shear modulus introduced phenomenologically 
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for solid quark matter with strangeness. More energy is then released, and 
this may enable a solid pulsar to glitch frequently with large amplitudes. Fur- 
thermore, the post-glitch behaviors may represent damped vibrations. In this 
paper, we try to model glitches in a starquake scenario of solid quark stars. 



2 The model 



A quake model for a star to be mostly sohd was generally discussed by Baym 
& Pines (flih and others who parameterized the dynamics for solid 

crusts, and possible solid cores, of neutron stars. Strain energy develops when a 
solid star spins down until a quake occurs when stellar stresses reach a critical 
value. Two possibilities were considered previously. Ruderman (till) assumed 
that the entire strain energy is relieved in the quake, while Baym & Pines 
fflih suggested only part of the stress is released and that the plastic flow is 
negligible. We present a third possibility that, during a quake, the entire stress 
is almost relieved at first when the quake cracks the star in pieces of small 
size (the total released energy may be converted into thermal energy -Etherm 
and kinematic energy E]^ of plastic flow, E^ = -Etherm + -^k) , but the part of 
might be re-stored by stress due to the anelastic flow (i.e., the kinetic energy 
is converted to strain energy again). As shown in Fig.l, a quark star may 
solidify with an initial oblateness Eq] stress (to be relative to Eq) increases as 
the star losses its rotation energy, until the star reaches an oblateness e+i. A 
quake occurs then, and the reference point of strain energy should be changed 
to El (the oblateness of a star without shear energy) at this moment. Damped 
vibrations in the potential field (see the inserted up-right part of Fig.l) occur 
therefore. After a glitch, the star solidifies and becomes elastic body again. 
No memory of the first glitch affects the second one. 

The density of quark stars with mass <~ l.SM© can be well approximated to 
be uniform (j2^. As a star, with an initial value Eq, slows down, the expected e 
decreases with increasing period. However, the rigidity of the solid star causes 
it to remain more oblate than it would be had it no resistance to shear. The 
strain energy in the star reads (18) 

-f^strain = B{e — Eq^ , (1) 



and the mean stress a in the star is 



1 dEstia 

V Oe 



fl{Eo - e) 



(2) 
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Fig. 1. The oblateness as a function of time, s{t). The horizontal sohd hnes are 
for the reference points of the strain. The steps, drawn very exaggeratedly, are 
starquakes. No starquake occurs anymore when e falls below (Tc//U (i.e., the stress 
can not be large enough for a quake). 

where e is stellar oblateness^ , V = AirR^/S is the volume of the star, and 
fi = 2B /V is the mean shear modulus of the star. We note that, in the following 
calculations, the stress-increase developed by the decrease of oblateness (due 
to the spindown) is only included. However, azimuthal stress due to the general 
relativistic effect ([llh (maybe similar to the frame-dragging effect in vacuum) 
of rotating solid stars may contribute significantly, though, unfortunately, the 
theoretical answers to elastic relativistic-stars with rotation are very difficult 
to be worked out. 

The total energy of such a rotating star with mass M and radius R is the sum 
of the gravitational energy -Egravi, the rotation energy Ej-ot, the strain energy 
-E'strain, the bulk energy E^, and the surface energy Eg, 

E = -£/gravi + -E'rot + -^'strain + E^ + E^ 

= Eo + As" + LV(2/) + B{e - Sif + E, + E,, (3) 

where Si is the reference oblateness before the {i + l)-th glitch occurs, Eq = 
—3M'^G/{5R), I is the moment of inertia, L = IQ is the stellar angular 
momentum, = 27r/P (P the rotation period), and the coefficients A and B 



^ The eccentricity e is defined by = 1 — c^/a^ (a and c are the semimajor 

and semiminor axes, respectively). The oblateness (or ellipticity) e = (/ — 

{Iq is the nonrotating, spherical moment of inertia) is related to e through e = 

(l_e2)-l/3_i^g2/3_ 
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measure the gravitational and strain energies (jisl ). respectively, 



A = (4) 

25 R ' ^ ^ 

B = ^rf/i. (5) 



The changes of and Es are much smaller than that of -Egravi, -E'rot, or -^strain 
when a star spins down, according to a mass formula for strange quark matter 
to be analogous to the Bethe-Weizsacher semi-empirical mass function (23). 
We therefore neglect E^ and Eg in the following calculations. 

By minimizing the total energy E, a real state satisfies (note that dl{e)/de = 
/o), 

Inn' B 



4(A + 5) A + B 

The reference oblateness is assumed, by setting i? = in Eq.(6), to be 

= hn?l{AA). (7) 

A star with oblateness of Eq.(7) is actually a Maclaurin sphere. When the star 
spins down to fi, the stress develops to 

, A hn? , 



'A + B " A{A + By' ' ' 

according to Eq.(2). A glitch takes place if the stress is greater than a critical 
one, a > ac- 

The first quake is characterized by a sudden shift of e, the amount of which is 

Ae = £+1 - e^i, (9) 

where the reference point is changed from Eq to Si. Due to the conservation of 
stellar angular momentum, the sudden change in the oblateness results in an 
increase of spin frequency, 

Afi A/ , 

- = -- = A.. (10) 

Starquakes would also produce a thermal energy dissipation during glitches 
that could be expected to be observable as an increase of X-ray luminosity 
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Fig. 2. The fractional increase in the spin frequency during a glitch, Ail/fi, vs. the 
rotation period, P. The lines labelled I, II, and III are for different shear modulus, 
fi. 1: fi = W^^ erg/cm^, II: fj, = 10^^ erg/cm^, and III: fj, = 10'^^ erg/cm^. The 
dash-dot, dotted, solid and dashed lines are for Uc = 10^^, cJc = 10^^, cJc = lO^'^ 
and o"c = 10^^ erg/cm^, respectively. We observe that the curves tend towards the 
asymptotic line when cJc <~ 10^^ erg/cm^. The stellar mass and mean density are 
chosen to be mass M = IAMq and density p = 4 x lO^'^ g/cm^ in the calculations. 

soon after glitch. According to the conservation of energy, one obtains from 
Eq.(3), 

O.^e-iml + B{e-i - ei)^ + ^thermal, (11) 

where Qi, which can be obtained by cr = cTc from Eq.(8), is the spin frequency 
when the first quake occurs, -E'thermai is the released energy in a starquake. 
The observed glitch size AQ/Q can be calculated from Eq.(6) to Eq.(ll). 
The calculations, based on these equations, are parameterized by five input 
quantities: M, p, /i, cXc and -^thermal, and two of them are fixed to be M = 
IAMq and p = 4 x 10^^ g/cm^. 

Fig. 2 shows the results for various /i and ac- We see that, for >~ 10^^ 
erg/cm^, AQ/Q decays to be a constant for pulsars with large rotation peri- 
ods. Small fi may generally result also in a small jump of AQ/Q for a certain 
period P. However, for cXc <~ 10^^ erg/cm^, the asymptotic line is reached, 
which means that ln{AQ/Q) oc — In P. We choose £'therm = 10^^ erg/cm^ in 
computation, but the curves do not depend sensitively on -Ethcrm- It is found 
that the glitch amplitude, AQ/Q, could be as high as observed if the shear 
modulus jj, ~ X0^°~'^^ erg/cm'^. Though the quark matter with low temper- 
ature and high baryon density is focused recent years, it is still impossible 
to determine the properties of such QCD phase by first principles (most of 
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the calculations are started from QCD phenomenological models). One of the 
proposed states of strange stars at low temperature could be a solid state (jl^ 
(and also the 2nd paragraph of §4), but the calculations on solid state is much 
more difficult than that of fluid one. Nevertheless, we may estimate the shear 
modulus n originated only by electric interaction between charged n— quark 
clusters {n: the quark number in a cluster) and a uniform background of elec- 
trons, through a similar quantum mechanical calculation of metals (01). The 
effective shear modulus, averaged over polarizations and directions, can be 
well fitted by /i ~ 0.12N{Ze)'^ /a in asymptotic case, where Z is the charge 
of quark-clusters, the cluster number density and a the separation between 



clusters ()25r ). For strange quark matter with baryon number density ub and 
electron number density to be ~ 10~^ that of quarks, one comes to 



where no = 0.16 baryons per fm^ is the nuclear saturation density. This result 
can be regarded as a low limit modulus of solid strange quark matter since 
the van der Waals-type color interaction, with a high coupling constant (to 
be corresponding to the electric charge e in electromagnetic interaction), may 
result in a larger shear modulus. We might then conclude that fi in the range 
of 10^'^"^'* erg/cm^ is not impossible. If the kilohertz quasi-periodic oscillations 
are relevant to the global oscillation behavior of solid quark stars, the shear 
modulus could be (|lCl ) fi ~ 10^^ erg/cm'^, which is much larger than the value 
(~ 10^^ erg/cm^) of neutron star crust. This could be unsurprise due to much 
high density. 

After the first quake, from Eq.(8), the stress in the star builds up again, at a 
rate of 

a = -f.e = -^^J^f^. (13) 



which is almost a constant for a certain pulsar with P and P during a period 
when the effect of P is negligible. Another quake occurs after a time of 

tq = (14) 



The calculation for certain P and P is presented in Fig. 3. For a certain set of 
{P, P}, the interval tq become longer for a higher value of a^. It is, however, not 
necessary to expect that a radio pulsar should jump its rotation periodically, 
because of no reason to show that the critical value cTc (and maybe /x) keeps a 
constant during its life. The time separation tq may present a quasi-periodic 
nature if cxc has a Gaussian distribution. 
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Fig. 3. Contour lines of the time separation between two quakes, tq, on the P — P 
diagram. This glitching interval time tq is labelled to the lines in unit of years. The 
parameters are taken as: M = \AMq, p = 4 x lO^^g/cm^, and /x = 10'^^ erg/cm^. 
Results with three values of o"c (10^^, 10^^, and 10^"^ erg/cm^) are for indications. The 
pulsar data are downloaded from http:/ / www . at nf . csiro.au/ research /pulsar /p sr cat ^ 

-10-' 



-10-= 

-10'* 
-10-^ 

-10-^ 



10'° 10'° 10" 10'° 

E.^ , (erq) 
thermal ' 

Fig. 4. The difference of e — i Si clS cl function of -Etherm- The lines are grouped for 
different periods; "F: P = 1 s, "11": P = 0.1 s, and "III": P = 0.01 s. The solid 
and dashed lines are for cjc = 10^^ and cjc = 10^^ erg/cm^, respectively, which are 
almost the same for certain period P. Other parameters chosen are: M = 1.4Mq, 
p = 4 X 10^^ g/cm^, and p = 10^^ erg/cm'^. 

3 The post-starquake behaviors 



II 



When i—th quake occurs in a pulsar, the oblateness decreases to at first, 
while the reference point is changed to Ei from ej-i. Soon after the rotational 
jump (e+j the star has a tendency to reach the equilibrium oblateness 

Ei. Certainly, this recovery process depends on the difference, E-i — Ei, which 
is shown in Fig. 4 as a function of i^therm- Increasing oblateness can obviously 
contribute an additional spindown role besides that of magnetodipole radia- 
tion, which results in larger P. Such a post-starquake could be observed as 
a recover behavior during post-glitch. From Fig. 4, one can see the values of 
E-i — Ei are generally not sensitive to both -Etherm and cTc. 

The recovery of E-i Ei is actually a complex process, in which both elas- 
tic revert and plastic flow could not be negligible. Nevertheless, this recovery 
process might be an analog of damped oscillations, with a mathematical de- 
scription of 

x + 2(3x + (jjIx = 0, (15) 
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with X = e — Ei. The second term in Eq.(15) arise from the anelastic effect 
during the recovery. The solution of Eq.(15) depends on the relative strength 
of the elastic spring and that of the damping. It is well known that three 
possible solutions exist; the so-called underdamped case: j3 < ujq, the critical 
case: (3 = ujq, and the overdamped case: /5 > c^o; where uq is the intrinsic 
frequency. One comes then to ~ \J jiR/M by dimensional analysis (force 
~ al^ and displacement ~ xR). 

In all the three cases, the oscillation damps by a factor of ~ exp(— with 
a typical time of r = 1/(3. We are just to discuss the critical damped motion 
below for an indication. In this case, c^o = 13 ■, the time evolution of e after the 
occurrence of starquake can be solved by Eq.(15), 

e{t) =ei + - ei){l + t/T)e-''\ (16) 



4 Conclusions and discussions 



A starquake model for pulsar glitches is developed in the regime of solid quark 
stars, and it is found that the general glitch behaviors (i.e., the glitch amplitude 
AVL/VL and the time interval tq) could be reproduced if solid quark matter 
has properties of shear modulus /x = 10^''~^^ erg/cm^ and critical stress cxc = 
]^gi8~24 QYg/cw? . It is suggcstcd that the post-glitch process could be described 
as damped oscillations, especially in the critical and the overdamped cases. 
Anyway, this is only a primary and simplified study of quakes in solid quark 
stars, more elaborate work, with possible modifications, on both quake and 
postquake processes is necessary in order to understand the nature of solid 
quark matter through glitching pulsars. 

We are dealing with solid quark stars in this paper. The quark Cooper pair- 
ing of the BCS type is suggested in quark matter of low-temperature but 
high baryon density, which may result in a color superconducting state (liel), 
with a large pairing gap on the order of 100 MeV. This kind of condensa- 
tion in momentum space takes place in case of same Fermi momenta; whereas 
"LOFF"-like state may occur if the Fermi momenta of two (or more) species 



of quarks are different (|27l). For 3 flavors of massless quarks, all nine quarks 
pair in a pattern which locks color and flavor symmetries, as called color- 
flavor locking (CFL) state ff28t) . However, for such quark matter, there exists 
a competition between color superconductivity and solidification, just like the 
case of laboratory low-temperature physics. One needs weak-interaction and 
low-mass in order to obtain a quantum fluid before solidification. This is why 
only helium, of all the elements, shows superfluid phenomenon though other 
noble elements have similar weak strength of interaction due to filled crusts 
of electrons. The strong color interaction (and the Coulomb interaction in 
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the system with strangeness) may be responsible for a possible solidification 
of dense quark matter with low temperatures. Further experiments (in low- 
energy heavy ion colliders) may answer whether quark matter is in a state 
of solid or color-superconductivity. Can a solid neutron star be possible? The 
answer might be no, because at least the part of neutron matter with approx- 
imate nuclear saturation density should be in a fluid state. In this sense, only 
solid quark matter is possible, and a quark star is identified if one convinces 
that a pulsar is in a solid state. 

We have assumed that the entire strain energy -Egtrain is relieved in a quake, 
which results in the reference oblateness of Eq.(7). However, it is possible that 
not entire, but most of, the energy -Estrain is released in a real situation, and 
the actual reference points are near but larger than £j. Of course, there is a 
tendency of £ — > £j after the i—th quake, but the effective shear modulus, /ies, 
of matter broken could be much smaller than that of perfect elastic solid, fi. 
The recovery timescale could be r ~ 15 days if fies is order of 10^^ erg/cm^. 

A large Vela glitch on 2000 January 16.319 was noted (|29l). and Chandra 
observations were earned out ~ 3^5 and ~ 35 days after tlTgl.tch Q. but 

no temperature change expected in conventional models with released thermal 
energy of ~ 10^^ ergs is detected. This could be understood in this starquake 
model since (1) the thermal conductivity of quark matter is much larger than 
that of hadron matter and (2) the thermal energy released, E'therm, to be much 
smaller than 10^^ ergs is possible. 
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